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Combining analytical and numerical methods, we investigate properties of the two-body ran- 
dom ensemble (TBRE) . We compare the TBRE with the Gaussian orthogonal ensemble of random 
matrices. Using the geometric properties of the nuclear shell model, we discuss the information 
content of nuclear spectra, and gain insight in the difficulties encountered when fitting the effective 
interaction. We exhibit the existence of correlations between spectral widths pertaining to different 
quantum numbers. Using these results, we deduce the preponderance of spin-zero ground states in 
the TBRE. We demonstrate the existence of correlations between spectra with different quantum 
numbers and/or in different nuclei. 
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I. INTRODUCTION 

The spectral fluctuation properties of complex nuclear 
spectra often agree with predictions of random-matrix 
theory, more precisely: With those of the Gaussian or- 
thogonal ensemble of random matrices (GOE). This is 
true for the resonances observed at neutron threshold 
and at the Coulomb barrier for protons but applies in 
a number of cases likcvifise to levels at lower excitation 
energies, see Rcf. 2] and references therein. 

The GOE is an ensemble of matrices where every state 
in Hilbert space interacts with every other state (if both 
carry the same conserved quantum numbers). Indeed, 
the probability density V{Hgoe) of the real and sym- 
metric Hamiltonian matrices 7?goe of the GOE has the 
form 

V{Hgoe) = AAexp { - ^TraceiHloE]} (1) 

where N denotes the dimension of the matrices and A 
a parameter which determines the average level density 
while A/" is a normalization factor. Eq. |^ shows that 
the matrix elements connecting any pair of states are 
uncorrelated Gaussian-distributed random variables. 

This basic tenet of the GOE is not in keeping with 
the shell model, the fundamental dynamical model of 
nuclear physics Indeed, that model is basically a 
single-particle model with a residual interaction. The 
interaction is dominated by two-body forces. In a repre- 
sentation where the many-body states are Slater deter- 
minants of single-particle states, a two-body interaction 
will have non-zero matrix elements only between those 
Slater determinants which differ by at most two units 
in the occupation numbers of the single-particle states. 
Out of the total number of such determinants, this is a 
small fraction. Put differently: In an arbitrary basis for 
the many-body states, the number of independent ma- 
trix elements of the two-body interaction is very much 
smaller than that of the GOE. This fact is changed only 
quantitatively but not qualitatively when we allow for a 



three-body residual interaction. 

Already in the I970's, this fact has led to the question: 
Are the predictions of the GOE for spectral fluctuation 
properties in keeping with the results of shell-model cal- 
culations with a residual two-body interaction? The an- 
swer ^ased on numerical calculations, has been affirma- 
tive [j, |5j , and numerous more recent calculations have 
confirmed it (see the review p and references therein). 
The calculations were based upon a random-matrix en- 
semble (the two-body random ensemble (TBRE)) which 
differs from the GOE and accounts for the specific prop- 
erties of the nuclear shell model: The existence of a resid- 
ual two-body interaction which conserves total spin, par- 
ity, and isospin. 

In spite of the agreement between GOE predictions and 
numerical results for the TBRE, open questions remain. 
In fact, analytical studies of the TBRE are virtually non- 
existent. Treating the TBRE analytically poses a severe 
challenge. We use a combination of numerical and ana- 
lytical techniques to display some of its properties. We 
focus on the underlying geometric structure of the TBRE 
and its implications for correlations between spectra of 
different quantum numbers (like nuclcon number and/or 
spin) . The results we are obtain are interesting also from 
a practical point of view, as we gain insight into the in- 
ner workings of the shell model and are able to separate 
geometric and dynamical properties. Together with the 
work reported in Refs. 7, 8], the present paper may be 
seen as a step towards filling the gap in understanding 
the TBRE. 

In treating the TBRE, wc neglect the influence of pos- 
sible three-body forces. It will become clear below that 
this omission affects our results in a quantitative but not 
in a qualitative fashion. Spectral fluctuation properties 
of the GOE type emerge only when there are at least 
three nucleons in a major shell. This is why wc confine 
ourselves to this case. 

This paper is organized as follows. The two-body ran- 
dom ensemble is defined and compared with the GOE in 
Section ^ Some of basic properties of the TBRE are 
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displayed in Section IlIII This main section contains re- 
sults about the information content of nuclear spectra, 
the preponderance of spin-0 ground states, and correla- 
tions between spectra with different quantum numbers. 
Section HVl contains a brief summary. Technical details 
are relegated to several Appendices. 



II. DEFINITION OF THE TBRE 

The two-body random ensemble (TBRE) is defined 
within the framework of the spherical nuclear shell 
model In that model, nuclcons move independently 
in a central potential with a strong spin-orbit force. We 
consider one of the major shells of that model. Our nu- 
merical examples are calculated for the sd-shell, with 
single-particle states labelled 51/2,^3/2: and and 
single-particle energies £1/2, £3/2, and £5/2- Our general 
considerations apply likewise, however, to other major 
shells in heavier nuclei. There, however, the number of 
many-particle states becomes forbiddingly large for nu- 
merical work. We sometimes also consider a single j-shell 
with half-integer single-particle total spin j. Although 
not realistic for nuclei, this case is sufficiently simple to 
yield useful insights. 

Putting several nucleons into a major shell, we con- 
struct a basis of orthonormal antisymmetrized many- 
body states of fixed total spin J, parity P, and isospin 
T. These states are labelled |J/j.), with J standing for the 
three quantum numbers J, P, T and with = 1, . . . , d(J) 
a running index with range given by the dimension c?(J) 
of Hilbert space 7i(J). We focus attention on a fixed but 
arbitrary z-projection M of total spin J so that d(J) is 
the actual number of states not counting their degener- 
acy regarding M . In the middle of the sd-shell and for 
low values of J, c?(J) is typically of order 10"^ and much 
larger for heavier nuclei (other major shells). The actual 
construction used for the basis of states |J/i) is immate- 
rial for what follows: The bases resulting from different 
modes of construction are connected by a unitary trans- 
formation. In the sd-shell, all single-particle states have 
positive parity, and it is not neccessary to carry the quan- 
tum number P. Likewise, we often consider only sd-shell 
states with isospin T = 0. In that case, it suffices to label 
the many-body states by the total spin J only. 

The number of nucleons in the major shell is denoted 
by m. Sometimes, we will consider simultaneously several 
nuclei with different values of to. In this case, we will 
denote, f.i., the dimension of Hilbert space by d{3,m), 
and similarly for other quantities. 

The many-body states |J/i) are eigenstates of the 
single-particle shell-model Hamiltonian with a very high 
degree of degeneracy. The degeneracy is lifted when we 
take account of the residual interaction of the shell model. 
We assume that this interaction mixes states only within 
the same major shell. That assumption is rather unreal- 
istic because intruder states from higher shells occur even 
at low excitation energies, and mixing with higher shells 



is bound to play a major role at the upper end of the spec- 
trum. The approximation is perfectly adequate, however, 
for the purposes of this paper. For the same reason we 
confine ourselves to a two-body residual interaction (al- 
though there is evidence that three-body forces do play a 
role in realistic attempts to model the spectrum) , and we 
omit the Coulomb interaction between protons. Finally, 
to focus attention on the role of the residual interaction, 
we assume that the single-particle energies within a ma- 
jor shell are all degenerate. For the sc?-shell, this means 
that we put £1/2 = £3/2 = £5/2 = 0. Then, the full 
Hamiltonian H of the shell model is determined entirely 
by the residual interaction. In the sequel we consider the 
matrix elements Hf^^^J) of that Hamiltonian taken with 
respect to the basis of states |J/i). 

Within a major shell, the residual two-body interac- 
tion V2 possesses a finite number of two-body matrix ele- 
ments. These have the following form. Let j;, i = 1, 2, 3, 4 
denote 4 (equal or different) values of total single-particle 
spin, parity, and isospin 1/2. Coupling ji and j2 (j3 and 
j4) to total two-body spin si (s2, repectively) and isospin 
ti (t2, respectively), denoting the parity of the resulting 
wave functions by tti {tt2, respectively), and introduc- 
ing the notation s for the quantum numbers s,7r, i, the 
reduced two-body matrix elements of V2 within a ma- 
jor shell have the form (j3j4s||V2||jij2s) where we have 
put Si = S2 = s since V2 conserves spin, parity, and 
isospin. The number a of such two-body matrix elements 
within a major shell is limited. For the sd-shell, we have 
a — 63, while for a single j-shell and identical nucleons, 
a = j -1-1/2. For brevity, we denote these matrix elements 
by Va with a = 1, . . . , a. We use the label a also for the 
specific two-body operator the two-body matrix element 
of which is Va, see Appendix 1. Within a major shell, V2 
is characterized completely by the a matrix elements Va, 
no matter how complicated the actual form of V2 ■ 

The Hamiltonian of the shell model is linear in the 
matrix elements Va and has the form 

H^,{J) = Y,Vo.C^,43,a) . (2) 

a 

The matrices C^^(J,a) transport the two-body interac- 
tion into the Hilbert space 7i(J) and depend upon the 
quantum numbers J, upon the particular states fj, and 
ly, and upon the particular two-body operator a con- 
sidered. The values of the C^,y(J,a)'s are completely 
specified by the underlying shell model, i.e., the single- 
particle states that occur within the given major shell, 
the coupling scheme used to construct the many-body 
states |J/i), and by the exclusion principle. The matrices 
Cfii, (J, a) do not depend upon the choice of the two-body 
interaction. That choice is specified by the values of the 
matrix elements Va- Eq. (j^J yields a decomposition of 
H into parts which are determined entirely by the sym- 
metries of the shell model (the matrices C^,y(J,a)), and 
parts which carry the information on the specific details 
of the two-body interaction (the matrix elements Va). 
We aim at generic statements about spectral properties 
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of H which apply to (almost) all two-body residual inter- 
actions. To this end we employ the TBRE: The matrix 
elements Va are assumed to be uncorrelated Gaussian- 
distributed random variables with mean values zero and 
a common second moment . Without loss of generality 
we put — 1 (we recall that all single-particle energies 
are equal so that the scale of the spectrum is determined 
by v^). Mean values of observables are worked out by 
integrating over the random variables Va, the measure 
being given by the product of the differentials of the Va 's 
and a Gaussian factor exp(— ?;^(a)/2). Having cal- 
culated the mean value and the square root of the vari- 
ance of an observable, we are sure that (within the error 
given by the latter) the mean value applies to all mem- 
bers of the ensemble, i.e., to all two-body interactions, 
with the exception of a set of measure zero. With the 
iiq's Gaussian random variables, the Hamiltonian _ff^^(J) 
represents an ensemble of Gaussian-distributed random 
matrices, the TBRE. 

As mentioned in the Introduction, numerical studies 
have shown that the spectral fluctuation properties of the 
TBRE generically coincide with those of the GOE. That 
result is tantamount to saying that H mixes the basis 
states |J/z) completely, irrespective of the specific choice 
of the iia's, and, thus, reflects a property of the matrices 
C^j/(J, a): Almost every linear combination of these ma- 
trices must be a sufficiently dense matrix in Hilbert space 
7i(J), with matrix elements which are sufficiently com- 
plex, to achieve such mixing. This statement is rather re- 
markable because the matrices C^y(J, a) are defined en- 
tirely in terms of an independent-particle model (which is 
intcgrable). In principle, the Cpy(J,Q;)'s can be worked 
out using group-theoretical methods! Intuitively, that 
mixing property of the C^y(J, Q!)'s can be understood by 
observing that each matrix element of C^,y(J, a) contains 
sums of products of Clebsch-Gordan and Racah coeffi- 
cients and coefficients of fractional parentage. For more 
than three particles in a major shell, the combination of 
these coefficients is highly complex, even though each one 
of them is well-defined and simple. A more detailed dis- 
cussion of the properties of the matrices C^^(J,a) was 
given in Ref. • Some further properties are displayed in 
Appendix 2. We show, in particular, that the C^^(J, a)'s 
do not commute. 

Before exploring some properties of the TBRE in the 
next Section, we compare the TBRE with the GOE. 
To this end, we recall some properties of the GOE, see 
Eq. (Q. The number of independent matrix elements of 
the GOE is N{N + l)/2. That is also the number of 
uncorrelated random variables which is, thus, very much 
larger than the dimension N of the matrices if ^ 1. 
Except for the symmetry of the GOE Hamiltonian iJcoE 
about the main diagonal, each matrix element carries an 
independent random variable. The situation is very dif- 
ferent for the TBRE. Here, the number a of independent 
random variables is typically much smaller than the di- 
mension (i(J) of Hilbert space. The complete mixing of 
the basis states and the ensuing validity of Wigner-Dyson 



statistics for the spectrum cannot be achieved by such a 
small number of random variables alone. In an essential 
way it is also due to the matrices Cp,y(J,Q!). These ma- 
trices form a fixed "scaffolding" for the TBRE. Although 
determined group-theoretically, these matrices are suffi- 
ciently dense and complex to guarantee complete mixing 
of Hilbert space for almost all choices of the random vari- 
ables Va- 

The GOE is invariant under orthogonal transforma- 
tions. Moreover, it is universal (i.e., the local spectral 
fluctuation properties do not depend upon the assump- 
tion that the ensemble has a Gaussian distribution) and 
ergodic (i.e., in the limit N ^ oo and for almost all 
members of the ensemble, the running average of an ob- 
servable taken over the spectrum of a single member of 
the ensemble is equal to the ensemble average of that ob- 
servable taken at fixed energy). These properties make 
the GOE to a formally attractive, mathematically acces- 
sible and universal model for stochastic spectral fluctua- 
tions. Does the TBRE possess any of these properties? 
The TBRE is manifestly not invariant under orthogo- 
nal transformations since the matrices C(J, a) are fixed 
by the underlying shell model. It is not clear whether 
the TBRE is universal. One might think, for instance, 
of modifying the TBRE in such a way that it incorpo- 
rates known features of the two-body interaction. One 
such feature is the pairing force. Doing so would im- 
pose certain fixed relations among the random variables 
Va- One could then assume that the remaining uncon- 
strained random variables have a Gaussian distribution. 
It is not known whether such a constrained TBRE would 
still generically possess spectral fluctuations of the GOE 
type. Finally, the ensemble is not ergodic. Indeed, er- 
godicity presupposes that the limit of infinite matrix di- 
mension can be taken. But the matrices C(J,a) of the 
TBRE have fixed dimension, and the question of ergod- 
icity cannot even be formulated meaningfully in this con- 
text. The only exception is the model of a single j-shell 
for which the limit j — s- oo is meaningful (it shares this 
property with French's embedded ensembles). It is not 
known whether in that limit, the TBRE for a single j- 
shell is ergodic. 

From all these points of view, the GOE is clearly the 
more attractive model for stochastic spectral fluctua- 
tions. Why then bother with the TBRE which is math- 
ematically so much more complicated and unrewarding? 
The TBRE offers the singular advantage of combining 
stochastic modelling with a realistic description of nu- 
clear spectral properties. This leads to speciflc insights 
which go beyond the realm of canonical random-matrix 
theory. Some of these are explored in the sequel. 



III. PROPERTIES OF THE TBRE 

For a full theoretical treatment of the TBRE, we would 
need a thorough analytical understanding of the proper- 
ties of the matrices C^,y(J, a). Unfortunately, we are still 
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far from this goal. We confine ourselves in this Section 
to a number of results which illuminate some properties 
of the TBRE. 



A. Information Content of Nuclear Spectra 

Almost by definition, a GOE spectrum carries no infor- 
mation beyond the symmetry of the Hamiltonian about 
the main diagonal. Indeed, the N{N + l)/2 matrix ele- 
ments of the Hamiltonian are drawn at random from a 
Gaussian probability distribution. Conversely, how many 
pieces of information are needed to reconstruct the GOE 
Hamiltonian from the spectrum? Knowledge of all the 
N eigenvalues does not suffice. Knowledge of one eigen- 
function adds — 1 additional pieces of information (all 
expansion coefficients with respect to a fixed basis except 
for the constraint due to normalization); knowledge of a 
second eigenfunction adds another N —2 such pieces (all 
expansion coefficients constrained by normalization and 
orthogonality to the first eigenfunction). Continuing that 
way, we see that the knowledge of all eigenvalues and of 
all eigenfunctions is needed to reconstruct the original 
GOE Hamiltonian. Since a GOE Hamiltonian does not 
carry any information, it is useless to investigate a GOE 
spectrum spectroscopically. 

Can we conclude from this result that whenever a spec- 
trum displays Wigner-Dyson statistics, it is useless to 
apply spectroscopy? In particular, would such a state- 
ment hold for nuclei? If not, which is the information 
content of a nuclear spectrum? These are the questions 
behind the title of the present Subsection. 

The answer to the first two questions is obviously a 
resounding NO. A Hamiltonian of the TBRE is specified 
by the choice of the a random variables Va- We recall 
that a is typically much smaller than d(J), the dimension 
of Hilbert space of the many-body states with quantum 
numbers J. Therefore, the knowledge of a ^ d(J) inde- 
pendent pieces of information from a spectrum of states 
with quantum numbers J alone will suffice to determine 
not only H{3) but also the Hamiltonian matrices H{3') 
of all other quantum numbers J' 7^ J! (This statement 
carries a small proviso: All the eigenvalues Sq(J) defined 
in Eq. Q below must differ from zero. Otherwise, addi- 
tional information from spectra with J' ^ J is needed). 
Given more than a pieces of data, the system of equations 
for the a matrix elements Va is obviously overdetermined, 
and nuclear spectroscopy is eminently meaningful even 
though the spectra may display Wigner-Dyson statistics. 
This conclusion is not confined to low excitation energies 
but applies, in principle, also to the domain of slow neu- 
tron resonances. The reason is that the mixing of basis 
states which leads to Wigner-Dyson statistics, is mainly 
due to the scaffolding matrices C^,y(J, a) which are fixed 
by the shell model. 

Can we identify the a independent pieces of informa- 
tion which are best suited to determine the variables Dq,? 
More importantly, do such a independent pieces of in- 



formation always exist? To answer these questions, it 
is useful to rewrite the TBRE in another form. We use 
the fact that such questions are most easily addressed in 
the framework of linear independence and orthogonality. 
The matrices C(J, a) span a linear space S. In S, the 
Hamiltonian |2Jl can be viewed as a vector with compo- 
nents Va- The space S can be endowed with a scalar 
product defined by the normalized trace. For two ma- 
trices C{J,a) and C(J,/3), that scalar product has the 
form 

SapiJ) = d-\3) Trace[C(J,a)C(J,/3)] . (3) 

The definition (jS)) has been used in a similar context 
before 0. With a, (3 — 1,2, ...,a, the scalar products 
Saf3{3) of all the C(J,a)'s form a matrix of dimension 
a. (We recall that in contradistinction, the C(J,Q!)'s 
themselves are matrices in the d(J)-dimensional Hilbert 
space H{3) spanned by the many-body states |J/i)). The 
C(J,q;)'s are real. The cyclic invariance of the trace 
then implies that the matrix Saf3{3) is real and symmet- 
ric. Moreover, the matrix Sap{3) is positive semidefinite. 
Therefore, Sa^iJ) can be diagonahzed by a real orthog- 
onal matrix Oap{3), 

{0{3)S{3)iOi3)fU = 5^,sli3). (4) 

The eigenvalues s„ (J) are non-negative and are arranged 
by decreasing magnitude, Si(J) > S2{3) > . . . s^{3) > 0. 
We define the matrices 



i3^,(J,a) = }_^0„^(J)C^,(J,/3) 
0=1 



and write Eqs. I0J and lO as 
1 



d(J) 



Trace[S(J,a)B(J,/3)] = Sa^p ^(J) 



(5) 



(6) 



For all s^(J) with s^(J) = 0, Eq. © implies that the 
trace norm of i?(J, a) vanishes which in turn implies that 
_B(J, a) = 0. In other words, there exist linear combina- 
tions of the matrices C{3,a) which vanish identically. We 
denote the number of non-vanishing eigenvalues s^(J) 
by ai, their positive square roots by Sa(J). The linear 
space S then has dimension ai. In iS, we define the oi 
orthonormal matrices 



5^1/ (J, a) = 



1 



Sa(3) ^ 



Oa/3(J)C^.(J,/3) for si(J)^0 



. (7) 

These are linear combinations of the matrices C^^(J,/3) 
and obey 

1 



d{3) 



Trace[B(J, a)B{3, (3)] = 5af3 with a, /3 = 1, 2, , 



, ai 
(8) 

The matrices B^^{3, a) form an orthonormal basis for the 
space S. Using these definitions, we rewrite the Hamil- 
tonian of the TBRE as 



H^A.3) 



E 

a=l 



Wa (J)Sa(J)5^^(J,Q;) 



(9) 
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The random variables Wa{3) are linear combinations of 
the Wq's, 

a 

WaiJ) =J20o.f}{3)Va . (10) 

/3=1 

The Wq(J)'s depend upon J via the orthogonal trans- 
formation 0{3). Like the v^s, the Wa(J)'s are uncor- 
related, have mean values zero and a common second 
moment equal to unity. The form © of the TBRE is 
totally equivalent to the original form ||2J). The number 
of variables Wa{J) appearing explicitly in the Hamilto- 
nian ^ is Oi and, thus, smaller than a if one or more 
of the eigenvalues s^(J) vanish. We shall see below that 
that is always the case. According to Eq. @ the trace 
norm of the Hamiltonian H{3) is given by 

^ Trace[i72(J)] = £ u;^(J)4(J) ■ (H) 

This expression is a measure of the spectral width of the 
spectrum of H{J). The spectral width is seen to be a 
random variable itself, with mean value ■s^(J). This 
last fact offers a physical interpretation of the eigenvalues 
s^(J) each of which gives the average spectral width of 
the spectrum due to the random variable z«q(J). 

Except for accidental degeneracies among the non- 
vanishing eigenvalues Sq(J), the matrices Bfj_^{J,a) are 
uniquely defined in terms of the diagonalization of 
'S'q/3(J). The orthonormality of these matrices is, of 
course, unchanged under orthogonal transformations of 
the linear space S, and it is possible to consider different 
forms of H arising from such transformations. Among 
these, the form is priviledged because the roles of 
the random variables Wa(J) and of the eigenvalues s'^{3) 
are well separated. In all other forms, these two entities 
become mixed by the orthogonal transformation. Never- 
theless, there exists another form of H which is physically 
interesting. This form is presented in Appendix 3. 

We return to the question raised at the beginning of 
this Subsection: Can we determine the variables of the 
Hamiltonian from a suitable set of data? With H 
given in the form of Eq. ((SJ, the relevant variables are 
now the Wq,(J)'s. Since the matrices B^y{3, a) have trace 
norm unity, the answer depends obviously on the eigen- 
values s^(J). If one or several of these quantities vanish, 
the Hamiltonian of the TBRE will, in fact, depend on 
a smaller number of random variables than its form (j^J 
suggests. And if one or several of the s^(J)'s are much 
smaller than the largest, then their influence upon the 
spectrum and eigenfunctions will be small, and it will be 
difficult or impossible to determine the values of the asso- 
ciated Wq's from a data set with experimental errors. We 
are, thus, led to study the distribution of the eigenval- 
ues Sq(J) or, equivalently, of their positive square roots 

Sa(J). 

We have calculated the square roots Sq(J) of the eigen- 
values s^(J) numerically for several cases: For all J~~ 



values for m = 6,m = 7 and m = 8 identical nu- 
cleons in the j — 19/2 shell, for all J- values of the 
T = states in ^^Ne, and for all J-values for the T = 
states in ^^Mg. We first address the case of the single 
j — 19/2 shell. Figures and |2] show the numerical 
results for the square roots Sa{J) for m = Q and for 
m = 8 identical nucleons, respectively. We observe that 
for all values of J, si(J) is distinctly larger than all other 
roots, and changes very slowly with J, increasing almost 
monotonically with increasing J. The same statement 
(near independence of J) applies to the corresponding 
eigenfunction Oia{J,m). We have calculated the over- 
lap function h{m,m'; J, J') = '^^Oia{J,m)Oia{J' ^m') 
both for m = m! and J ^ J' and for J = J' but 
m ^ m! . Results for the first case were shown in Fig. 4 
of Ref . 8] . These demonstrate that the overlap function 
decreases very slowly with increasing distance | J — J'|. 
Even for | J — J'| = 20, the overlap function is still larger 
than 0.8. As for the second case, we have calculated 
the overlap function for J = J' = 0,1,2, ...,10 and 
(TO,m') = (6,8), for J = 0, 1, 2, . . . , 10 with J' ^ J+l/2 
and {m,m') — (6,7), and for J' — 0,1, 2,..., 10 with 
J = J' -1- 1/2 and (m, m') = (7, 8). The overlap function 
is larger than 0.995 in all cases except for pairs involv- 
ing m = 6, J = 1. This last fact is due to a special 
feature of the matrix Saf3 for m = 6, J = 1: Two eigen- 
values of Sap are zero. We conclude that the matrices 
B^i,{J, TO, a = 1) consist of (nearly) the same linear com- 
binations of the matrices C^i/( J, m, a) for many values of 
J and TO. 



I I I 




10 20 30 40 



J 

FIG. 1: (Color online) The square roots of the eigenvalues 
of the matrix Sa/siJ) for a single j-shell with j = 19/2 and 
m — 6 identical nucleons versus total spin J. 

Some of these features can be understood analytically 
via arguments presented in Appendix 4. We use the nu- 
merically established near independence of the largest 
eigenvalue si(J) and of the associated eigenfunction 
Oia(J) on J to study the average Xap of SafS (averaged 
over all J~values). We do so in the limit (2j -1-1) ^ 1. 
In this limit, Xap factorizes, Xaf3 ~ f{ci)f{P)- Hence all 
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FIG. 2: (Color online)Same as Fig.Qbut for m = 8 nucleons. 

eigenvalues but one vanish, and is given by /^(a). 
For m ~ 6 {m — 8) this yields si 6.3 (si « 11.4, respec- 
tively). We also determine the eigenfunction Oia which 
(except for normalization) has components (1,5,9,...). 
The resulting matrix Bi is 

Bi cx^(2s + l)C(J,l + s/2) . (12) 

The sum runs over s = 0, 2, 4, . . . , 2j — 1 for the values of 
the two-particle spin. The matrix Bi is the matrix of the 
monopole operator Bi . Inspection reveals that this oper- 
ator is approximately proportional to the unit operator. 
Thus, all matrices that are orthogonal (under the trace) 
to Bi must have approximately zero trace. Moreover, the 
leading term in the transformed Hamiltonian ^ is ap- 
proximately a scalar. It defines the approximate value of 
the centroid of the eigenvalues and has presumably only 
a small influence on the spectral statistics. The value of 
the centroid is determined by the values of si(J) and of 
k;(1). The spectral statistics are largely determined by 
the remaining terms in H. 

Our analytical results are mirrored in the numerical 
findings. Except for the largest J-values, the matrices 
B(J, 1) are close to the unit matrix in d{J) dimensions. 
Comparison with Figs. and |21 shows that the analytical 
results for the eigenvalues are semiquantitatively correct. 
For the eigenfunction Oia the numerical diagonalization 
for m = 6 and J = yields an overlap of 0.97 with the 
analytical result. 

We turn to the data for the sd-shell. Our calculations 
are partly based on the shell-model code Oxbash 10] . 
The square roots Sq( J) for ^"Ne are shown in Fig.O those 
for ^*Mg are shown in Fig. 0] The tendencies are quali- 
tatively similar to those in the single j-shell. In ^"Ne, 
si(J) is not larger than all other roots by as big a factor 
as in the single j-shell, or as in ^"^Mg. This suggests that 
the clear separation of si(J) from the rest of the roots 
requires a minimum number of nucleons. Indeed, in Ap- 
pendix 4, the factorization of Xaf3 occurs only for m > 4, 



and Si grows with m roughly like m^. Thus, our results 
suggest a generic behavior for many nucleons in the shell. 
Again, we find that the leading matrix B{J, 1) is (except 
for a normalization factor) approximately equal to the 
unit matrix. 
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J 

FIG. 3: (Color online) The square roots of the eigenvalues of 
the matrix Sai3{J) for the sd-shell nucleus '^"Ne with isospin 
T = versus total spin J. 
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FIG. 4: (Color online) Same as Fig.^lbut for ^''Mg. 

Another remarkable feature in all the figures is the 
disappearance of one or several roots for some values of 
J. In fact, for each spin, at least one eigenvalue vanishes 
identically. The corresponding scalar two-body operator 
is given by — J{J+ 1) where J denotes the operator of 
total spin. The matrix representation of this operators is 
a linear combination of the matrices C^^{3). In the case 
of the sd-shell nuclei, an analogous statement applies to 
the isospin operator T^ — r(r-|-l). More eigenvalues than 
one vanish in cases where the number of independent 
matrix elements of H{3) is smaller than the number of 
two-body matrix elements. 
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How do these results relate to the actual determination 
of the effective shell model interaction from experimental 
data? We recall that usually one proceeds in two steps. 
In a first step, the effective interaction is determined 
through a G-matrix calculation that is based on nucleon- 
nucleon potentials which fit nucleon-nucleon scattering 
data with very high accuracy. The resulting G-matrix 
usually does not reproduce experimental nuclear spectra 
well enough, and corrections have to be made in a second 
step. There are basically two ways to make these correc- 
tions. A first approach insists on minimal corrections and 
only corrects the monopole term by a fit to experimental 
data. This is very reasonable since the monopole term 
has the largest trace norm, and is essentially the only 
matrix with a finite trace (or centroid). This last prop- 
erty of the monopole has been known since the pioneering 
work of Pasquini and Zuker IjJ . We note that the result- 
ing shell-model Hamiltonian has an impressive predictive 
power, see, e.g., the shell-model studies of /p-shell nuclei 
by the Strasbourg group [l2| . In the second and alterna- 
tive approach, one attempts to fit all two-body matrix 
elements to experimental data. This is done through an 
iterative procedure that starts from the G-matrix and 
makes small changes in the matrix elements until a best 
fit is obtained. As the number of data points (i.e., energy 
levels) usually exceeds the number of two-body matrix 
elements, the problem is over-determined, and only a 
smaller number of linear combinations of data points can 
be used in the fit due to rank-deficiency problems of the 
matrix equations involved ^13j . As an example we men- 
tion the widely used Brown-Wildenthal interaction 
for sd-shell nuclei where only 47 linear combinations of 
the 66 parameters (63 two-body matrix elements and 
three single-particle energies) are determined by fit while 
the remaing ones are taken from the G-matrix. In larger 
model spaces the ratio of well-determined linear combi- 
nations usually decreases. For the /p-shell, only 70 linear 
combinations out of 199 free parameters are used in the 
fit to data As shown by our results in this Sub- 

section, the answer to the question "Which are the most 
important linear combinations of two-body matrix ele- 
ments" is less determined by the data used in the fit but 
rather by the inherent structure of the shell model itself. 
Our work identifies which linear combinations these are: 
The Wq(J)'s which correspond (for fixed J) to the largest 
eigenvalues. Unfortunately, the property of an eigenvalue 
to be large or small changes little (if at all) with J, and 
the same is true of the orthogonal transformation defin- 
ing the random variables WaiJ)- This fact shows that 
even an enlarged data set might not necessarily lead to 
more accurate values for the u;ct(J)'s. 

In conclusion, the information content of a nuclear 
spectrum as given by the TBRE is very different from 
that of a GOE. The actual possibility to determine the 
unknown parameters Wq-(J) depends on the magnitudes 
of the eigenvalues Sa(J)- We stress that these eigen- 
values are uniquely determined by the underlying shell 
model and do not depend upon the residual interaction. 



Thus, they may be calculated prior to any attempt to 
fit the residual interaction to actual data. We believe 
that such a procedure might be useful for practical shell- 
model work. 



B. Preponderance of Ground States with Spin Zero 

In 1998, Johnson, Bcrtsch, and Dean found that 
the TBRE is likely to yield spin-zero ground states in 
even-even nuclei with a probability which is consider- 
ably larger than the fraction of states with spin zero in 
the total shell-model space. This came as a surprise be- 
cause the TBRE does not have a built-in pairing force 
or quadrupole force. Subsequent work showed that sim- 
ilar regularities exist in bosonic ^3 and electronic [l^ 
many-body systems with two-body interactions. The 
phenomenon of spin-zero preponderance seems a robust 
and rather generic feature which has received much atten- 
tion since, see the reviews 0, and references therein. 
Here we discuss it in the framework of our representa- 
tion lO of the Hamiltonian of the TBRE. We extend our 
earlier work in Ref. 8]. As before in this paper, we fo- 
cus attention on the case of a single j-shell and on the 
T — states of nuclei in the sd-shell and, therefore, use 
the label J rather than J. 

Our approach is based upon the following considera- 
tion. For each value of J and each realization of the 
random variables we define the spectral radius R{J) 
as the distance of the farthest eigenvalue from the center 
of the spectrum. (This can be the smallest or the largest 
eigenvalue. In view of the randomness of the signs of 
the Wa{J)^s, that ambiguity is irrelevant). The R{J)'s 
are random variables. For every value of J we determine 
the probability that the spectral radius R{ J) is maximal. 
To this end we relate R{J) to the spectral width cr(J). 
The square of the spectral width a{J) is given by the 
normalized trace norm of the Hamiltonian, see Eq. Hll|l . 

'^'(.J) = f2^liJ>liJ) ■ (13) 

a=l 

We postulate a linear relationship between the spectral 
radius and the spectral width, 

R{J) = r{J)a{J) . (14) 

We show that the scaling factors r{J) are nearly con- 
stant (i.e., independent of the random variables Wa{J)) 
but do depend upon J. Then, Eq. (|14|l relates the ran- 
dom variables R{J) with the random variables a{J). The 
dependence of the scaling factors r( J) on J and the dis- 
tribution of and correlations among the spectral widths 
cr(J) together lead to an understanding of the prepon- 
derance of spin-zero ground states. 

We first address the scaling factors r{J). The inset of 
Fig-Elshows the dependence of R{0) on (t(0) for 6 identi- 
cal nucleons in the j = 19/2 shell, together with a linear 
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fit, for 900 realizations of the TBRE. We see that Eq. ^ 
holds with an approximately constant value of r(0). The 
reduced per degree of freedom is about 0.9. This value 
increases toward 1.5 for spins around J « 20. Figure |S1 
shows the linear-fit values of r( J) versus J. These have 
an overall tendency to decrease with increasing J. This 
fact reflects the overall decrease of the number d{J) of 
states of spin J with increasing J. Indeed, the smaller 
the number of states in the spectrum the closer do we ex- 
pect i?( J) and (t( J) to be. This expectation can be quan- 
tified. We recall that a shell-model spectrum with spin J 
calculated from the TBRE has approximately Gaussian 
shape ■ For a Gaussian spectrum with spectral width 
a normalized to the total number d of states, the average 
level density p{E) has the form 



piE) 



2tt(t 



(15) 



The most likely position of the lowest (the highest) state 
in the spectrum is found by integrating p{E) from —00 
(H-cxd, respectively) to the point where the area under the 
integral equals 1/2. This point is equal to R. Thus, 



/27r 



dx exp(— a;^/2) 



(16) 



The solutions R/cr — r{d) of this equation depend on 
d. In Figure 13 we have displayed the resulting values for 
r{d{J)), again for the j = 19/2 shell with m = 6 identical 
nucleons. The odd-even staggering reflects correspond- 
ing changes in the dimensions d{ J). The error in our 
theoretical determination is expeced to be of the order of 
the mean level spacing, i.e., of order d~^ and, thus, small 
for d 3> 1. The actual discrepancy between the data and 
our analysis is much larger. We ascribe this to the fact 
that the assumed Gaussian shape of the spectrum holds 
only approximately. 

We turn to the distributions of and the correlations 
among the spectral widths cr(J). From Eq. (|13|l . the 
ensemble average (henceforth denoted by an overbar) of 



KJ) is 



and the variance is 



a4(J)- (a2(J))' = 2^4(7) 

^ ^ a=l 



(17) 



(18) 



The normalized root of the variance equals \j2ja\ when 
all eigenvalues are equal and equals \pl when all eigen- 
values but one vanish. Therefore, we conclude from Fig- 
ures ^ to 0] that the fluctuations of cr^(J) are biggest for 
the largest J-values. This is also borne out when we 
calculate the complete probability distribution Pj{cr) = 
2a6{a'^ — o'^(J)) for (t(J). We write the (5-function as a 




FIG. 5: (Color online) Dependence of the scaling factor r(J) 
in Eq. 1141 on spin J for m — 6 identical nucleons in the j = 
19/2 shell. Solid line: Data based on the linear fit. Dashed 
line: Prediction based on Eq. 11611 . Inset: Dependence of i?(0) 
on cr(0) for 900 random realizations of the TBRE. Solid line: 
The linear fit which determines r(0). 



Fourier integral and perform the integrations over all the 
Gaussian variables Wa{J). We find 



dt e' 



^arctan(2tsQ, (J)) 



n(l-f 4^2,4 (J))l/4 



(19) 



Plots of Pj{u) for various values of J were shown in Fig- 
ure 3 of Ref. 8]. The fiattest curve is the one for the 
largest value of J as expected. 

Central to an understanding of the preponderance of 
spin-zero ground states in nuclei are the correlations be- 
tween the spectral widths u{J) pertaining to different 
values of J. We have seen that both the eigenfunc- 
tion pertaining to the largest eigenvalue and all eigen- 
values of the matrix Safi change little with J. In the 
extreme case where for a pair { J, J') of spins we would 
have s^(J) = s„( J') and Wa{J) = Wa{J') for all values of 
a, the spectral widths would obviously be totally corre- 
lated: Their joint probability distribution is proportional 
to a delta function S{a{J) — a{J')). In this case, the value 
of the level with the lowest spin depends not on a{J) but 
only on r(J). Actually the eigenvalues are not exactly 
equal and, more importantly, the eigenfunctions for the 
smaller eigenvalues differ, and so do the corresponding 
i/7c((J)'s. Still, we must expect significant correlations. 

The correlations cannot be worked out directly from 
Eq. ifT^ as this would require knowledge of the correla- 
tions among the ( J)'s for different values of J. Rather, 
we use that ct^ is defined in terms of the trace norm of 
the Hamiltonian, use Eq. ((2Jl for the latter and the defi- 
nition for Saf3{J) and have 



a/3 



(20) 
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For the covariance of two a s this yields 



(72(J)a2(J')-CT2(J) ct2(J') = 2 Trace[S'(J)S'(J')] ■ (21) 

In general, the right-hand side of this equation differs 
from zero. It is not possible straightforwardly to calcu- 
late the joint probability distribution of cr(J) and o'(J')- 
However, the probability 



p(J,J') = e(a2(J)-a2(J')) 



(22) 



that spin J has larger spectral width than spin J' can be 
calculated in terms of the eigenvalues qa of the matrix 
{Sap{J) — SafsiJ'))- The calculation is quite similar to 
the one which leads to Eq. H19|l and yields 



p( J, J') 



1 



dt sin(i J2i3 arctan2%) 



(23) 



This formula is in good agreement with our numerical 
results. 

We emphasize that our approximations for r{J) and 
(t( J) do not yield a faithful representation of the distri- 
butions and correlations of the spectral radii R{J). How- 
ever, these approximations turn out to be sufficiently ac- 
curate to predict the probability that R{J) = r{J)a{J) 
has maximum value. These predictions are compared 
with the results of numerical diagonalization in Figs. 
toO In all four Figures, we plot the probability that the 
ground state has spin J versus J. The data points are 
calculated from 900 realizations of the TBRE. The solid 
lines show the probability that the spectral width cr(J) 
has the largest value. The dashed lines show the proba- 
bility that the product r{J)a{J) has the largest value. 




FIG. 6: (Color online) The probability that the ground state 
has spin J versus J for m = 6 identical nucleons in the j — 
19/2 shell. Dots: Data points from 900 realizations. Solid 
line: Probability that the spectral width o"( J) has the largest 
value. Dashed line: Probability that the product r{J)a{J) 
has the largest value. 

We observe that the states with spin zero and with 
maximum spin have the largest spectral widths. The 




FIG. 7: (Color online) Same as Fig.^lbut for m = 8 identical 
nucleons 
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FIG. 8: (Color online) Same as Fig.|Slbut for the T = states 
in 2"Ne. 



factor r(J) suppresses the high J-values. It causes a 
staggering in the probabilities which is not present for 
the spectral widths. In all Figures, the agreement of the 
dashed lines with the data points is satisfactory. We note 
that the predictions are somewhat less accurate for ^''Mg. 
A more detailed analysis of this case shows that the de- 
termination of the scale factor r(J) by fit is somewhat 
less accurate as the typical P^r datum is around 2.5. 

It is of some interest to see whether our arguments 
apply also to nuclei with an odd number of nucleons and 
half-integer total spin. In complete analogy to Figs. 
to |51 we have calculated the probabilities for the half- 
integer spin states for m = 7 identical nucleons in the 
single j — 19/2 shell, to form the ground state. The 
results are shown in Fig. 1101 with symbols that carry 
the same meaning as before. Again, the agreement is 
very convincing. This strengthens our belief that our 
approach is generic. 

In summary, we have shown that our semi-analytical 
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FIG. 9: (Color online) Same as Fig.|Slbut for the T = states 
in ^-^Mg. 




FIG. 10: (Color online) Same as Fig.|S|but for m = 7 identical 
nucleons 



approach to the TBRE offers a satisfactory explanation 
for the preponderance of ground states with zero spin 
observed for this ensemble. We have established an ap- 
proximate proportionality between two sets of random 
variables, the spectral radii i?( J) and the spectral widths 
cr(J), R{J) ~ r{J)a{J), with nearly constant scaling fac- 
tors r{J). We have presented semi-analytical results for 
r(J), and closed expressions for the distribution func- 
tions for the spectral widths as well as for some of their 
correlations and joint distribution functions. The scaling 
factors r{J) depend essentially on the dimensions d{J) of 
the associated Hilbert spaces, show odd-even staggering, 
and an overall decrease with increasing J. The spectral 
widths reflect properties of the underlying shell model 
and fluctuate most strongly for large values of J. Their 
average values are largest for J = and for maximum 
J. In the products R{J), the large J-values are sup- 
pressed by the scaling factors r(J). The good agreement 
of our theoretical predictions with the data in all cases 



suggests that we have identified the generic mechanism 
which is responsible for the preponderance of spin-zero 
ground states. That mechanism should likewise operate 
in other major shells and perhaps in other many-body 
systems governed by two-body interactions. 



C. Correlations between Spectra with Different 
Quantum Numbers and/or in Different Nuclei 

For different quantum numbers J ^ 3' and/or different 
nucleon numbers m ^ m', the Hamiltonian matrices of 
the TBRE defined in Eq. ||2Jl depend upon the same set 
of random variables and are, therefore, correlated. This 
fact is rather obvious from the point of view of nuclear 
physics: Changing the residual interaction will simulta- 
neously change all the spectra in all the nuclei within 
the major shell under consideration. However, from the 
point of view of random-matrix theory, the existence 
of correlations between spectra each of which displays 
Wigner-Dyson type level fluctuations, is excluded by as- 
sumption. To the best of our knowledge, a theoretical 
framework for the treatment of such correlations does 
not exist. The parametric level correlations which have 
been discussed earlier, in the context of both condensed- 
matter theory 22] and random-matrix theory do 
not cover the present case. Indeed, there one considers a 
Hamiltonian which depends upon an external parameter 
like the strength of an external magnetic field. The di- 
mension of the Hamiltonian matrix remains unchanged as 
the parameter is varied. Here, in contradistinction, cor- 
relations exist between Hamiltonian matrices which differ 
in matrix dimension. Moreover, parametric level corre- 
lations generically tend to zero as the difference between 
the old and the new values of the parameter increases. In 
the present case, it is not clear whether the correlation 
between two spectra differing, say, by two units of total 
spin is larger or smaller than that between two spectra 
differing by ten units. Also, it is not clear whether such 
correlations are big enough to be significant experimen- 
tally. These facts prompt us to explore the magnitude of 
spectral correlations both versus spin and versus mass. 
We focus attention on the case of sd-shell nuclei. 

Using the standard deflnition for the level density 



p{E,J)=J2S{E-E^{J)) 



(24) 



in terms of the eigenvalues Ef^(J) pertaining to spin J, 
we have calculated numerically the density-density cor- 
relations for ^"^Mg for T — states with J — and ,7=1. 
Our ensemble average uses 400 realizations of the TBRE. 
The results are shown in Fig.^] The left panel shows as 
a density plot the mean value of the product of the level 
densities, p{E, 0)p{W, 1), versus the energies E and W in 
units of The midddle pan el show s a nalogou sly the 
product of the mean densities p{E,0) x p{W,l). This 
is essentially the product of two Gaussian distributions. 
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The rig ht panel shows t he actu al correla tor, i.e., the dif- 
ference p{E, 0)p{W, l) — p{E, 0) X p(W, 1). The maximum 
value of the correlator is about 10% of the product of the 
mean level densities. This demonstrates the existence of 
weak correlations. The correlator has a minimum in the 
center of the two spectra. We cannot exclude the possibil- 
ity that with increasing matrix dimension, this minimum 
widens, so that significant correlations would exist only 
in the tails of the spectra. We have found correlations 



between spectra pertaining to different J-values also in 
the case of a single j-shell, where they are significantly 
stronger and are of the order of 50%. Correlations are 
also expected to exist between states in different nuclei 
when these are governed by the same TBRE. A case in 
point are the T — spin zero states in ^°Ne and ^"'Mg. 
Their correlations are shown in Figure IT^ where we use 
the same display as in Figure 1111 Again, we find weak 
correlations of the order of several percent. 




FIG. 11: Correlations between the J = 0,r = states and the J = 1,T — states of ^'*Mg (from 400 realizations of the 
ensemble). Left panel: Average of the product of the two level densities. Middle panel: Product of the averages of the two 
level densities. Right panel: The correlator. 




FIG. 12: Same as Fig.IHlbut for the J = 0,T = states of ^""Mg and of ^^Ne (from 400 reahzations of the ensemble). 



In view of the results in Subsection IIII Al and in Ap- 
pendix 3, the existence of density-density correlations is 
expected. Indeed, the form of the TBRE Hamilto- 
nian shows that the centroids of the spectra must show 
particularly strong correlations. It is essentially these 
correlations that are displayed in Figs. ITTI and 1121 We 
ask: Are there also significant correlations between the 
spectral fluctuations of levels with different spins and/or 
different masses? These would show up in correlations 
between level spacings and would, thus, be independent 
of the fluctuations of the centroids of the spectra. To 
answer this question, a finer test than displayed in Figs. 
111! and [T^ is required. We have not explored this ques- 



tion any further. We recall, however, that data relevant 
to this question were pubhshed in Ref. O: With Ej 
the energy of the lowest state with spin J, the ratio 
p ~ {E4 — E2)/{E2 — Eq) shows a broad peak in the 
range p — to p ~ 1. This suggests the existence of level 
correlations which go beyond those due to the centroids. 

As remarked above, spectral fluctuations of spectra 
with different quantum numbers go beyond the assump- 
tions of canonical random-matrix theory. Therefore, it 
would be of interest to verify the existence of such corre- 
lations experimentally. We have explored the possibility 
of such a test in the following way. In nuclei the residual 
interaction is fixed and not random. Therefore, ensemble 
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averages cannot be taken. In canonical random-matrix 
theory, one circumvents this problem using ergodicity: 
The running average of an observable over a single real- 
ization of a spectrum is equal to its ensemble average. In 
the TBRE, we propose to take the average over an ensem- 
ble of nuclei in the same shell as a substitute. In the spirit 
of this proposal (which unfortunately lacks a theoretical 
foundation) we have calculated the spacings between the 
lowest levels for a number of sd-shell nuclei using the 
Brown- Wildenthal two-body interaction The en- 

semble consists of the nuclei ^°~^''Ne, ^^"^"'Na, ^"'"^^Mg, 



26 Al, 30Si, 34p, 32,34g^ g^^^ 36 J^j. p^j. ^j^g ^^^^ (g^j^) ^^^^ 

nuclei of this ensemble we considered the correlations of 
the J = and J = 2 (J = 1/2 and J = 5/2) states. We 
label the nearest-neighbor spacings of levels with equal 
spins consecutively by AEi, i = 1,2,3, starting from 
the ground state. Figure [T51 shows . from left to right, the 
average of the spacings AEiAEj, the product of the av- 
erages AEi AEj, and the correlator AEiAEj — AEi AEj, 
in miits of (MeV)^. We see that correlations exist that 
are about 10% of the product of the averages. 
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FIG. 13: Level spacing correlations for 17 sd-shell nuclei between spins J = and J — 2 {J = 1/2 and J — 5/2 )for even-A 
(odd- A) nuclei. Left: Average of product of spacings. Middle: Product of average spacings; Right: Correlation. 



The results of this subsection show that small corre- 
lations (of the order of 10%) exist between the spectral 
fluctuations of sd-shell nuclei. We expect that the exper- 
imental verification of the existence of these correlations 
is somewhat challenging, as it requires the measurement 
of several complete spectral sequences involving 5-10 lev- 
els each. 



IV. SUMMARY 

We have studied various aspects of the two-body ran- 
dom ensemble in nuclei, and presented three main re- 
sults. First, the geometric properties of the nuclear shell 
model become transparent once one replaces the usual 
two-body operators by linear combinations that are or- 
thogonal under the trace. This transformation yields 
one operator that is approximately proportional to the 
unit matrix and has the dominant spectral width. This 
monopole operator sets the scale for nuclear binding, and 
can be derived analytically. The remaining linear combi- 
nations of two-body operators have approximately zero 
centroid and smaller spectral widths. They determine 
the spectral fluctuation properties. Some of these oper- 
ators have very small or even zero widths and therefore 
are close to or equal to the null operator. It is difhcult 
or impossible to determine the corresponding linear com- 



binations of two-body matrix elements by fit to experi- 
mental data. Thus, the geometric properties of the shell 
model explain the difficulties encountered when residual 
interactions are determined by fit. 

Second, we presented analytical and numerical results 
regarding the fluctuations of and correlations between the 
spectral widths of shell-model operators. The total spec- 
tral width can be linked to the spectral radius through a 
simple scale factor. This approach allowed us to give a 
semi-quantitative explanantion for the preponderance of 
spin-zero ground states both for a single j-shell and for 
the sd-shell. 

Third, we studied correlations between spectral fluctu- 
ations of spectra belonging to different spins and/or mass 
numbers. Our numerical results for the sd-shell indicate 
that correlations exist at the 10% level in the region of 
low-lying excitations. This is new territory and goes be- 
yond the assumptions of canonical random-matrix the- 
ory. It would be of considerable interest to explore how 
these correlations depend on matrix dimension and, thus, 
change as we consider other major shells. 

We have shown that the complex geometric structure 
of the shell model can be understood in rather simple 
terms. Apart from their theoretical interest, our results 
might be of practical use in fitting effective interactions, 
and they may motivate an experimental test of the cor- 
relations between spectra belonging to different quantum 



13 



numbers. 



Appendix 1. Two— body Interaction 

For simplicity of notation, we confine ourselves to the 
case of a single j~shell and identical nucleons, omitting 
isospin and parity quantum numbers. The generalization 
is straightforward. 

Let at and a,, be the creation and annihilation op- 
erators for a particle with half-integer total spin j and 
z-component /i. Two such identical particles are coupled 
to total spin s with s — 0,2, . . . ,2j — 1, 

= ^' " '")%aA/-M • (25) 



This is a tensor of rank s which is non-zero if s is even. It 
transforms under rotations with the complex conjugate 
Wigner I?-function, 



(26) 



N 



The adjoint operator is 

This operator transforms with the usual L'-function, 

{All)new = Y^{Dljpj){Alf)old ■ (28) 

N 

From the unitarity of the /^-functions, it follows that 

T{s) = ^(Al,)tAi, (29) 



M 



is a scalar. 

It is convenient to rearrange the order of the creation 
and annihilation operators, 

T{s) = - ^ c{jjs; ^, Ad - fJ.)c{jjs; v,M -v) 

fj.Mi' 



xatf ,,a,,a\.a 



2j + l^ t 



(30) 



With m the number of particles in the shell, the last term 
is equal to m{2j + l)/(2s + 1). This term is denoted by 



A{s) and gives a diagonal contribution to the Hamilto- 
nian. We define the tensor operators y^* of rank t 



Ail = E '^(■^■?^! ~ p)«p"p-ji/ 



(31) 



For t ~ 0, the tensor operator A*^ is just the number 
operator, for i = 1, it is the operator of total spin, etc. 
Using standard transformations of the Clebsch-Gordan 
coefficients, we rewrite the terms in Eq. H26II as 



r(s) = (2. + 1) ^ iy(jjjj; St) Y,i-r+'^+^'A\,A'_ 



M 



M 



-A{s) 



(32) 



With these definitions and with a — 1 + s/2, the two- 
body interaction V2 has the form 



j+1/2 
V2 = E ^<^^i^) 



(33) 



a = l 



where we use either the definition (|29|l or the defini- 
tion (|23) for T{s). 



Appendix 2. Some properties of the matrices C^„(J) 

As in Appendix 1, we confine ourselves to the case of 
a single j-shcU and identical nucleons, omitting isospin 
and parity quantum numbers. We use the definition H25|l 
and 1(23 . 

To work out the properties of the matrices C^j/(J), 
we shall presently need the commutator [T{s),T{t)] for 
s ^ t, with T(s) defined in Eq. (|29|l . That commutator 
is given by 



M.N 



[T{s),T{t)] = (iAl,)HAl„{A%)M% 



+ (A%)H{Al,)\A%]A%] .(34) 



The commutator of two A's can be expressed in terms of 
the irreducible tensor operators A introduced in Eq. I|31|l. 
Wc find 



[iAl,)\A%] = -4^V(2s + l)(2p+l)iy(,sjtj;jp) 

p 

xc{stp;M,N)Blj_j, . (35) 
Here is an irreducible tensor of rank p defined as 

B^n = E(^)^~"'^(-?-?P' ^ ~ '")4a„_™ • (36) 

n 

Hence, 
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cistp; M,N){A'j,)^Bl,^^Al, 



M.N 



■E 



(2f+l)(2p+l) 



W{tjsj;jp) 



2s + 1 

Y,c{tsp-N,M){Alj)'^Bl_,,A\ 

(37) 



M • 



M.N 



The result is a scalar because the three operators A'^ , i? 
and A are coupled to rank zero. It is obvious that the 
commutator \T{s),T{t)\ differs from zero and does not 
lie in the linear space spanned by the operators T. A 
recoupling of the operators shows, in fact, that the com- 
mutator is a sum of three-body operators. 

The elements of the matrices C^i/( J, s) with s = 2a — 2, 
may be viewed as matrix elements of operators C{J,s). 
The latter are defined in the space of Slater determinants 
of single-particle shell-model states (i.e., states not cou- 
pled to total spin J) and have the form 



C{J,s)^V{J)T{s)V{J) 



(38) 



where 'P{J) is the orthonormal projector onto the sub- 
space of many-body states with fixed total spin J. The 
projectors obey 

V{.h)V{J2) - V{.h)5j,j, ; [V{J)]^ = V{J) . (39) 

With the help of the operator J of total spin, V{J) is 
explicitly written as 



v{j) = n 



j{j +1)- 



(40) 



The multiplication extends over all possible values Ji of 
total spin. 

It is straightforward to show that the operators 7^ (J) 
commute with T{s). Therefore, we have 

C( J, s) = V{J)T{s)V{ J) = T{s)V{J) = V{J)T{s) . 

(41) 

It follows that the commutator [C{J,s),C{J,t)] can be 
written as 



[C{J,s),C{J,t)]=V{J)[T{s),T{t)] 



(42) 



But we have shown in Eq. (|37|l that the commutator 
[r(s), T{t)\ does not vanish and is, in fact, a sum of three- 
body operators. It then follows from Eq. (|42() that the 
operators C do not commute. Rather, their commuta- 
tors represent sums of three-body operators (projected, 
of course, onto the space with fixed J). In particular, it 
is not possible to diagonalize the Cs simultaneously. 



For a = 1 , . . . , a we define 

c(J, a) = Trace C^^(3, a) 
and introduce the traceless matrices 



(43) 



^^.(J, a) = C7^,(J, a) - ^^^y^-^;.. • (44) 

With these definitions, the Hamiltonian in Eq. ((2Jl takes 
the form 



i?M-(J)=iy^.(J) + ( E 



c(J, a) 
-^(Jf' 



Q 5t,u (45) 



where 



(46) 



is traceless. We can now diagonalize the real symmetric 
positive semidefinite a-dimensional matrix 

Sa.p = ^Trace [C(J, a)C(J, (3)] (47) 

by an orthogonal transformation Oapi^). Of the result- 
ing eigenvalues s^(J, a), exactly (ai — 1) differ from zero. 
Using the same construction as in Eqs. Q to ||SJ and 
denoting the resulting oi — 1 orthonormal matrices by 
i3(J, a), we arrive at 

ai — 1 

i?p.(J) = E '^'a(J)5a(J)B(J,a) 



+ -^M-E 



c(J, a) 

i^(jr' 



(48) 



Similarly as in Eq. (Uni but with O replaced by O, the 
(fli — 1) random variables Wc[(J) are linear combinations 
of the fa's. They are uncorrelated, have mean values 
zero, and a common second moment equal to unity. The 
matrices i3(J, a) are traceless and, hence, orthogonal to 
the unit matrix in the sense of the trace norm. Thus, 
Eq. H48|) is another representation of the TBRE in terms 
of oi orthonormal matrices. In contrast to Eq. lO, we 
have now separated a term (proportional to the unit ma- 
trix) which defines the centroid of the spectrum but has 
no bearing on the actual level sequence, from the re- 
mainder which in turn possesses a vanishing centroid but 
yields the actual details of the spectrum. The form H48|) 
is probably not useful for data analyses but should be 
helpful whenever we are interested in spectral fluctuation 
properties. These depend only upon the ai — 1 variables 



Appendix 3. Another Form of the TBRE 



Appendix 4. Approximate DiagonaHzation of Sa 



As mentioned in Subsection IIII Al there exists another 
possibly useful form of the TBRE. This we now derive. 



As in previous Appendices, we perform the calculation 
for a single j-shell with identical nucleons. The number 
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of nucleons is denoted by to. We use the fact that here 
we have s = 2a — 2 and replace a, /? by s, t. We use the 
definition H41|l and the fact that P{J) and T{s) commute 
to write Sst in the form 



SstiJ) 



1 



d(J) 



Trace [ViJ)T{s)T{t)] 



(49) 



The calculation is difficult for fixed total spin J. How- 
ever, our numerical calculations show that the eigenval- 
ues of SstiJ) ^^'^ the eigenvector pertaining to the largest 
eigenvalue change little with J for the lowest values of J. 
Moreover, these lowest values of J have the largest val- 
ues of d{J). Therefore, we approximate the calculation 
of Sst{J) by averaging over all values of J, thus defining 



X, 



1 



EjdiJ) 



(50) 



We observe that the sum over J removes the projectors 
V{J) since their sum equals unity. Thus, the average 
is equal to the normalized trace of T{sT(t) taken with 
respect to all states in the Hilbert space spanned by nor- 
malized Slater determinants (without any attention paid 
to the coupling of spins) . We denote that trace by angu- 
lar brackets and have 



Xst - {Tis)T{t)) 



(51) 



We calculate Xst approximately by keeping only the lead- 
ing terms in an expansion in inverse powers of (2j + 1). 
That is the meaning of the approximate sign in this Ap- 
pendix. 

The trace in Eq. (|51|l can be worked out using Wick 
contraction. The relevant term is 



N- 



(52) 



Two types of Wick contractions occur: (a) Those where 
contractions connect pairs of operators separately within 
the group of the first four and of the second four oper- 
ators, and (b) those where at least two operators of the 
first group of four are contracted with the same number 
of operators of the second group. Among all these, the 
terms of leading order in an expansion in inverse powers 
of (2j -|- 1) are the terms with the maximum number of 
unconstrained summations over magnetic quantum num- 
bers. Every Wick contraction of a pair of operators pro- 
duces a Kronecker delta. However, the two Wick contrac- 
tions in case (a) aff'ecting the first (the second) group of 
four operators produce the same Kronecker delta's (5^^ 
{Spa-, respectively). This is not so in case (b) where each 
of the two Kronecker delta's imposes a new constraint. 
We conclude that the terms of leading order arise only 
from Wick contractions for case (a). This fact implies 
that 

Xst « (T(s)) {T{t)) = f{s)f{t) (53) 
where we have defined 

f{s) = (T(.)) . (54) 



Factorization of X as in Eq. H53|) implies that the leading 
eigenvalue sf defined in Eq. Q is given by 



(55) 



while all other eigenvalues are w 0. The associated ma- 
trix Bi is the matrix representation of an operator Bi. 
The latter is proportional to X^s /(^)-^(^)- Normaliza- 
tion requires {Bf) — 1. We use again that the Wick 
contractions approximately factorize and find that 



EsPis) 



(56) 



We observe that (Bi) w 1. The same statement holds, 
within our factorization approximation, for all higher 
powers of Bi. Thus, Bi is (except for normalization) 
approximately equal to the unit matrix. With the di- 
mension of total Hilbert space given by N{j, m) = : 
this gives 



B^^ N-'^\j,m) 1 



(57) 



It remains to work out the function /(s). We denote 
the Wick contraction by an index W. We have 



t t 



Thus, 



[5fj.v — SpM-v_ 



W 



t t 

O-fj.O-M-fj.'^n'^M-ti ■ 

(58) 



f{s) - (T(s)) « E '^(■?'-?'*; M - Mjs; v,M-v) 

2 



/(«) 



2j + l 



2j + l 



2 (2s -f 1) . 



(59) 



(60) 



For j — 19/2 and m = 6 (to = 8), this yields for sf the 
value 39.933 (126.21) so that si = 6.32 (11.24, respec- 
tively) . 

With /(s) given by Eq. (|60|) . we can work out Bi ex- 
plicitly using Eq. (|56|l and Eq. (|30|l for T{s), omitting 
normalization factors. This yields 



Bi cx 



1 



-(2j + 1)2to1 



(E.(5i+l)2)i/n2^ 
- E (2s + l)c(jjs; M - Ai)c(jjs; M - u) 



SM^LV 

t t 



(61) 



We retrieve the unit matrix as in Eq. H57|) and con- 
clude that compared to the laeding term, the last sum 
in Eq. (EU is of order l/(2j + 1). 
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These results allow us to estimate sf (J). We do so by 
assuming that the results obtained by averaging over all 
values of J apply separately to each J. Since B{J, 1) is, 
except for terms of order l/(2j + 1), equal to a multiple 
of the unit matrix, we equate the variance of the coeffi- 
cient multiplying B{J, 1) in Eq. with the variance of 
the coefficient multiplying d^'^/'^{J)5^i, in Eq. We 
obtain 

s\{J)^{l/d{J))Y,c\j,a) . (62) 

a 

By the same token, we expect that all the ai — 1 eigen- 
values s'^{J) in Eq. 147|l are on average of equal size and 
of order l/(2j -I- 1). This implies that the statistical fluc- 



tuations of the centroid of the spectrum are considerably 
larger than those of individual level spacings. 

In summary, we have shown that the leading eigenvalue 
and eigenfunction of Xst are easily obtained because for 
2j + 1 ^ 1, Xst approximately factorizes. The factors 
f{s) and f{t) are simply given by Wick-contracting the 
operators T{s) and T{t) separately. The leading eigen- 
value is given by Eq. H55(l and the matrix Bi is approxi- 
mately given by the unit matrix. 
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